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Constraints on the Combined Models with ^ Inflation 
and Viable f{R) Dark Energy 

Chao-Qiang Geng^’^’^, Chung-Chi Lee^, and Shan 
Lin3 


Abstract We investigate the observational constraints 
on the modified gravity, which combines the in¬ 

flation with the power-law (exponential) type of the 
viable f{R) dark energy models. We discuss the dif¬ 
ference between the combined model and grav¬ 

ity in the inflationary epoch and obtain the constraints 
on the deviation power q as well as the parameters in 
f{R) by using the CosmoMC package. The allowed 
ranges of the spectral index and tensor-to-scalar ratio 
from the Planck data are highly restricted, resulting in 
q < 2.66 X 10“^ and 2.17 x 10“^ for the power-law and 
exponential types of f{R) gravity, respectively. 
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1 Introduction 


Various cosmologi cal observations such as Type-la Su¬ 


perno vae (SNIa) l|Riess et al.l Il998t iPerlmutter et al 


unive rse to provide a grac eful solution to these prob¬ 
lems ( Lvth fc Riottol 1998 ). As a result, our universe 
has gone through two different eras of the accelerating 
expansion: inflation and late-time acceleration. 

It is well known that the simplest way to realize 
the late-time acceleration is to introduce a cosmological 
constant into the Einstein equation. The correspond¬ 
ing theory is referred to as the ACDM model, which 
fits with the observational data very well, and becomes 
the standard model of cosmology. However, there exists 
several unnatural properties in this model, such as coin¬ 
cident and hierarchy problems. In order to remove these 
problems, people have tried various methods. One of 
which is to introduce a homogeneous and isotropic en¬ 
ergy density with a negative pressure into the theory of 
General Relativity (GR), and t he other one is to mod - 
ify the Einstein’s gravity theory ( Gooeland et al ][2^ . 
Among them, f{R) gravity is one of not only the min- 


199S ), cosm i c microwave 


__ imal exten sion of the AGDM model, but also the most 

background radi ation (GMB) ( Sper^f-Miftk^ eory to explain inflation and dark energy. 


2003 . 20061 : iKomatsu et al.l I2OO8I 2010) and b aryor 


acoustic oscillation fBAOll Eisenstein et al. [ioosf) indi¬ 
cate that our universe is currently under an accelerating 
expansion, which is known as the dark energy prob¬ 
lem. On the other hand, there are some cosmological 
problems in which the hot Big Bang model cannot be 
explained, such as the horizon, flatness and monopole 
problems. A rapid accelerating expansion era, known as 
inflation, is introduced right after the beginning of the 
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To describe the late-time accelerating phenomenon, 
several viable f{R) viable models have been proposed to 
satisfy the constraints from the oretical considerations 
and c osmological observa tions (jPeFelice fc Tsuiikawa 


2 OI 0 I: Bamba et al. 2011 1. These viable f{R) mod¬ 


els ca n be classified i nto power-law an d exponential 


2012 : Geng et al. 2014[l . Under these 


types (|Lee et al 
viable conditions, the future singularities have been 
removed and the instability of cosmological pertur¬ 
bation never occurs in the universe. Additionally, 
there still exists a broad window of experi mental tests 
by s pace gravitational wave detections (jYang et al 


[mH). Through the Markov chain Monte Garlo method 
(MGMG) analysis, the viable f{R) dark energy models 
get a smaller likelihood than that in the AGDM model, 
denotin g that they are preferred by cosmological obser- 


vations (Motohashi et al 


2m.3l:lGeug ellblbmdll . Be¬ 


sides the dark energy problem, /(i?) gravity is also used 
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to address the dark matter p ro blem dCembranos 2009 . 
2011 : Sharif &: Yousaf 2014al fbl. 2015c] bi as well as the 


2005fl . Since the MGC AMB program (Mo dification 


inflation scenario. In this manuscript, we will focus on 
the viable f{R) models with dark energy and inflation. 

It is well-know n that the inflat ion, first proposed 
by Starobinsky ( Starobinskv 198^. successfully de¬ 


scrib e s the reheat ing period (IVilenkinlll985t iMiiic et al. 


I986I : Ford I986ll and predicts a tiny tensor-to-scalar 
ratio, which is consi stent with the observations to- 
Since both f{R) gravity and 


2015). 


day (|Ade et al. 

R^ inflation achieve their goals by modifying the action 
of general relativity to a function of Ricci scalar, it is 
tempting to consider their actions together to provide a 
unified framework for both inflation and the late-time 
acceleration epoch. In fact, the R^ model remains one 
of the most attractive models today describing inflation 
because it solves the singularity proble ms when com¬ 


bined with f(R) dar k energy models (lAppleby et al 
20 I0HLee etS]l2012h . 


As a simple general i zation of t he R^ inflation, many 


people (iMartin et al. 

Rinaldi_^^_alj20Mc b d: Motohashi20I4 : 

2015: Noiiri fc Odintsovl200,1 . 


2014aH3: ICodello et al. 2014: 


2006 


2007 


of Growth with GAME) ( Hoiiati et al.l[2011 1 uses the 
parametrized framework to include f{R) gravity into 
GAME, we only consider the linear perturbation and 
assume that the background evolution is the same as 
the AGDM model. 

This paper is organized as follows. In Sec. O we first 
give a brief review on the R^~'^ inflation with slow-roll 
parameters, and then introduce the f{R) modifications 
in the scalar and tensor perturbations. We also esti¬ 
mate the spectral index rig and tensor-to-scalar ratio 
r. In Sec. [3l we show our results of the observational 
constraints on the combined models by using the Gos- 
moMC package. Finally, we present our conclusions in 
Sec. m 


2 The Action and Parameters 

Let us begin by considering the 4-dimentional general 
/(i?) action 


M^Tzakul et^J^ 


2008 


Gognola et ^ 


i\/f2 r 

^ J d'^x^/^f{R) + Sm , 


2008; Ehzalde^^alJ 201T ;jBamba et aL 20081 


( 1 ) 


2014; Eamba 


fc Odintsovl2015 


Ghakravartv fc Mohantd 120141) have investigated the 
possibility of the power in R^ to deviate slightly from 
two, named the R?~’^ inflation, where q <C 1. The R^~'^ 


inflation has been shown to be compatible with obser¬ 
vations and predicts a larger scalar-to-tensor ratio than 
that in the R^ inflation. It is thus natural to combine 
the R?~^ inflation model with viable f{R) dark energy 
models to ensure the late-time acceleration of our uni¬ 
verse. As an example, a recent attempt of combining 
these two has been given in Ref. ( Artvmowski fc LalakI 
201411 . In this paper, we would like to study the con- 


Ghakravartv et al.||2013t[ where Mpi = (87rG)“^/^ is the reduced planck mass, 

f{R) is a general function of the Ricci scalar i?, and 
Sm is the matter action, including relativistic and non- 
relativistic components. In this study, we consider f{R) 
to be the form, 


f{R) = R + F{R) + 


R^-<i 

M2-29 


( 2 ) 


straints on the combined inflation with dark energy 
models. The action includes the early universe 
fiati on and late-time f {R) dark energy of the St arobin¬ 


sky ( Starobinskvll2007tl and exponential gravity (IZhana 
2005 : Gognola et al.l l20o5 Ihinder 2009 : Eamba et al 


2010 |) models, which belong to the power-law and ex- 


where ^nd F{R) are responsible for the inflation 

and late time dark energy, respectively. It is worth to 
mention that the R? model was claimed to not only 
reali ze inflation but a l so sol ve the dark matter prob¬ 
lem ( Gembrano j 2009 . 2011 1. The later has also been 
discussed with the extended dark m atter m odels in the 
literature ( Sharif fc Yousaf l2014alli3.[20T5lbl) . 

In gener al, the FfR) models, which sa tisfy the viable 
conditions ( DeFelice fc Tsuiikawa 2010l) can be cate¬ 
gorized into two classes, the power law and exponen¬ 
tial types. To investigate the generic features in these 
classes, we concentrate on t he exponential type with the 


ponential types, respectively. 

We modify the Gode for A nisotropies in the Mi¬ 
crowave Background (GAME) ( Lewis eL jlJll999) and p-ypnnpTitial f orm o f gravity ( Zhan dl2005 : Gognola et al. 
the C osmological MonteCarlo (CosmoMC) package ([Lewis fecBst^deil 120091: iBamba et al.lboinll an d the power- 
20021) to study the constraints on inflation as well law o ne, called the Starobinsky model ( Starobinskv 
as f{R) dark energy models from the cosmological 
observat ions, including those of the C MB data from 


2007D. 


Planc k (|Ade et al 

201 2 ll 


201,111 and WMAP ([Hinshaw et al. 


BAG data fro m Baryon Oscillatio n Spectro- starobinsky : F{R) = -\R, 


scopic Survey (BOSS) ( Anderson et al.l[2012l) and SNIa 
data from Supernova Legacy Survey (SNLS) ( Astier et al.^ 


1 - 


1 


Exponential : F{R) = —/3i?c(l — , 


(3) 

(4) 
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where A and /3 are the model parameters and Rc rep¬ 
resents the constant characteristic curvature. 

2.1 inflation 




Fig. 1 The potentials in the Einstein frame as functions of 
the scalar field 4> in the combined model, f{R) = R+F{R) + 
R^/M^, where F{R) has the form of (a) the Starobinsky 
model with (A, n) = (2, 2) (solid line) and (4, 2) (dash- 
dotted line) and (b) the exponential gravity with (5 = 2 
(solid line) and /? = 4 (dash-dotted line), while the dashed 
line corresponds to the R^ inflation, f{R) = R + R^/M^. 


First of all, let us focus on the inflationary stage. By 
conducting a conformal transformation = fngfii, 

and fa = df{R)/dR = transform the 

/(i?) action from the Jordan frame into Einstein frame, 


5' = 


/ 




~ 1 

-!/(</.) 


(5) 


where R is the curvature in the Einstein frame and the 
potential is given as 


VW 


Mgi Rfn - f 


2 


fk 


( 6 ) 


Note that in the combined model, the potential depends 
on not only the term but also the F{R) dark 

energy model, and the scalar field in the Einstein frame 
is defined by 


to derive the explicit form of V{4>) in the combined 
model. In Fig. [T1 we numerically reverse Eq. d?]) to ob¬ 
tain the relation R = i?(^) and show the potentials as 
functions of the scalar field by substituting R{(j)) into 
Eq. ©with the R^ inflation and combined model. Note 
that, as an illustration, we simplify the model with 
(g, M'^/Rc) = (0, 10"^) and F{R) is the Starobinsky 
model in Eq. ([3]) (left panel) with (A, n) = (2, 2) (solid 
line) and (4, 2) (dash-dotted line), and the exponential 
gravity in Eq. (right panel) with (3 = 2 (solid line) 
and 4 (dash-dotted line).Obviously, when the curvature 
is large enough and close to the inflationary stage, the 
potentials of the combined models (solid line) mimic 
that of the R^ inflation (dashed line). In the small 
curvature regime (^ < 0+), the potentials of the R^ 
inflation and the combined models increase to infinity 
simultaneously. In the slow-roll regime, the R^~‘^ infla¬ 
tion model is no difference from the combined model. 
However, the post inflationary history in the combined 
models could have a significant change of that in the R^ 
inflation. On the other hand, in the post inflationary 
epoch, the curvature oscillation in the Jordan frame 
also exhibits a different behavior between the R^ in¬ 
flation and the combined models. In the R^ inflation, 
the scalaron oscillates around the potential minimum 
at ()> = 0 in the Einstein frame, corresponding to the 
curvature R oscillating around i? = 0 in the Jordan 
frame. On the contrary, the curvature oscillation in the 
combined models “never” reaches zero and evolves to 
be negative. As shown in Fig. [21 one has 

V {(t>{R))\R^Rcr- ^ ^ (8) 

where (A, n) = (2, 2) for the Starobinsky model 
(dashed line) and (3 = 2 for the exponential one (solid 
line).We choose the exponential F{R) gravity to be the 
example. The potential in the Einstein frame is given 
from Eq. ®. When the curvature is small enough and 
i?/M^ becomes negligible, V{(()) approaches infinity at 
Jr — 3 + Ffi —>■ 0, i.e., 

R —>• Rcr = Rc In (3 , (9) 

which represents a positive critical curvature. 


/* = eV^«"- = i + (2_,)^ + A. (7) 

From Eqs. ©, o and 0, we see that the large cur¬ 
vature regime corresponds to a large 4>, so that the 
behavior of the scalar field is dominated by the R^~‘^ 
term, whereas the small curvature regime represents 
a small or negative 0, resulting in the F{R) term re¬ 
sponsible for the evolution of the scalar. Due to the 
non-linear form of Eqs. and o, it is not easy 


As discussed above, the post inflationary epoch his¬ 
tory might be changed by including F{R) dark energy 
component. However, matter and dark energy compo¬ 
nents, corresponding to Sm and F{R) in Eq. Q still 
can be ignored during the inflationary stage, respec¬ 
tively, so that the inflation potential is given by substi¬ 
tuting f{R) = R + aR'^~^ into Eq. ([B]), 

V{(()) = K)e Mpi _ I 


( 10 ) 
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Fig. 2 The relation of the potential in the Einstein frame 
and the curvature in the Jordan frame, where the solid and 
dashed lines correspond to the combined models with expo¬ 
nential {j3 = 2) and Starobinsky (A = 2 and n = 2) gravity, 
respectively. 


in both Einstein and Jordan frames are the same, i.e. 
Ne — N. Substituting Eq. (ITTl) into Eq. m. we find 


q^O: 

N{cl>) 


q = 0: 
Nicjy) 


3(^ln 


_l_ 2(1 — g) 

n ^end 

qgV 3 Mpi _|_ 2(1 _ q'j 



( IT*— fT*—\ 

I ^ "pi — eS ^ "pi j 

f2cj)end\ 

4 3 Mpi Vs Mpi j ■ 


(19) 


( 20 ) 


where Vb = (1 — <?)(2 — q)^^~'^'^^^Mpi/2 is a 
mass dimension-four constant. 

Given the potential V{<j)), we obtain the slow-roll 
parameters, defined by 


e = 





With these parameters, we can write down the inflation 
observables, such as the scalar and tensor spectral index 
(ris^nt), tensor-to-scalar ratio (r) and scalar spectral 
index running (ag = dus/d\nk) to be 


Ug — 1 = —6e -1- 2ri , 

(12) 

nt = -2e, 

(13) 

r = 16e, 

(14) 

ag = 16e77 — 24e^ — 2^ , 

(15) 

respectively. The end of inflation is 
satisfy the condition, 

usually assumed to 

— 4^end ^ 

(16) 


Therefore, one obtains 

^end _ j \i]_ 

Mp)" “ V 2 

Now, we can derive the number of e-foldings during 
inflation. 


(2 + V3)(l-g) 
V3-(1 + V3)q 


N = 


r 


Hdt ~ 


Ml 


I 

J cb 


PI 4^end 


V 


d^, 


(18) 


where the tilde denotes the variable in the Jordan 
frame. Erom the conformal transformation, we have 
Hdi = Hdt 1 -I- ~ Hdt during the slow- 

roll regime. As a result, the numbers of the e-folding 


In order to avoid the negative value inside the logarithm 
in Eq. cu, we concentrate on g > 0 in this study. We 
can reverse the equations to get 

^ = ^iq,N). (21) 


The slow-roll parameters are given by 


e 




[qfR + 2{l-q)]^ 

Hl-q)HfR-ir' ^ ’ 

2[gVI-(l-g)(2-5g)/H + 4(l-g)^] 

3(l-g)2(/fl-l)2 ’ ^ ^ 

4[qfR + 2{l-q)][q^f% 

+ (l-g)(l-2g)(2-5g)/| 
-il-q)\lO-17q)fR 

+8{l-q)^][9{l-qf{fR-l)^]-\ (24) 


where Jr = ^ numerically ana¬ 
lyze the inflation by using Eqs. (fTOl) - (E^ . In 

Fig. [3^, we depict the scalar spectral index Ug (solid 
line) and tensor-to-scalar ratio r (dashed line) as func¬ 
tions of the power deviation q in the R^~’^ inflation 
model, where the black and gray lines represent the in¬ 
flation e-foldings with = 50 (black) and 70 (gray), 
respectively. In Fig.|3lD, we plot (j) as functions of q with 
N = 50 and 70 to demonstrate the values of the scalar 
field at the initial time as functions of q by Eq. (l3T|) .In 
Fig. [5J:, we show r as a function of rig, where where 
the long-dashed, solid and dashed lines denote N = 50, 
60 and 70, while the contour plot presents the ACDM 
model constrained at 68% and 95% confidence levels 
with the same datasets used in Sec. [3l respectively. 
Comparing these two plots, we observe that the R^~‘^ 
inflation allows a larger value of r, but q < 0(10“^). 
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Fig. 3 ( a) The scalar spectral index Us (solid line) and 

tensor-to-scalar ratio r (dashed line) as functions of the 
inflation power parameter q in the R^~'^ inflation model, 
where the black and gray lines correspond to = 50 and 
70, respectively, (b) The scalar at initial time, cj>, as func¬ 
tions of q with = 50 (black line) and 70 (gray line), (c) 
The tensor-to-scalar ratio r as a function of the scalar spec¬ 
tral index Us, where the long-dashed, solid and dashed lines 
denote N = 50, 60 and 70, while the contour plot presents 
the Icr and 2a bound in the ACDM model, respectively. 


2.2 F{R) Dark Energy 


After inflation, the scalaron settles down in the global 
minimum of the potential. Hence, the term in the 
general action is lost its importance and can be ignored 
in the evolution afterwards. We then have the same 
action as the F{R) action, which is 

r 

S=^ J d^x^{R + F{R)) + Sm- (25) 

One can obtain the field equation by varying the action 
with respect to the metric, given by 


FrR^u - ■^F{R)gij,i, + (g^i,n - V^iVv)Fr 




(26) 


where Fr = dF{R)/dR, □ = is the 

d’Alembertian operator and is the energy-momentum| 
tensor defined by 


T, 


fliy 


-2 SS^ 


(27) 


To deal with the linear perturbation in /(i?) grav¬ 
ity, let us first write down the Friedmann-Robertson- 
Walker (FRW) metric in the Newtonian gauge, 

= a^(T) [—(1 -I- 2T')dT^ + (1 ~ 2^)6ijdx’'dx^] , (28) 

where a(r) is the scale factor, r is the conformal time, 
and fk and are the scalar perturbations. The per¬ 
turbed energy-momentum tensor is 

T° = -ipM + SpM), (29) 

= —{pM + PM)vM,i , (30) 


where vm is the veloc ity held. By follow ing the simi¬ 
lar procedure in Ref. l|Hoiiati et al.N2011l) . we can de¬ 
rive the relationship between the metric potential to the 
perturbation of the matter density in the sub-horizon 
limit. 


a2 

$ 

¥ 


T' = -47rGe//(fc, a)pMSM , 


= l{k,a) , 


(31) 

(32) 


where Sm = ^pm/pm + SFIavM is the gauge-invariant 
matter density perturbation and 


r< ^ 4- qfiiAKH 
Geffik,a) = — - 

Fr 


(33) 

(34) 


We incorporate the linear perturbation effects by in¬ 
troducing the effective gravitational constant Ggj / and 
the ratio of Newtonian gauge scalar potentials as de¬ 
fined above. 

For completeness, we consider not only the scalar 
Newt onian potentials but also the tensor perturba¬ 
tion ( Zhou fc Hel 2014 ). The perturbation of the FRW 
metric is given by. 


ds^ = a^ij) [—-|- (dy -f Dij)dx''dx^^^ , (35) 


where Dij is a traceless divergence free tensor held. In¬ 
serting Eq. (IM|) into the held equation (l26ll . one can 
deduce the scale independent modified equation for the 
tensor mode evolution. 


D'l 


2)H + ) D', + = 


2 T 

mI,Fr ’ 


(36) 


where R = a~^da/dT, the prime denotes the derivative 
with respect to the conformal time, and ir'k is the tensor 
perturbation of Ty. 
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3 Constraints from cosmological observations 

In this section, we use the current observational data 
to constrain the inflation as well as the model 

parameters in F (R) with dark energy. We perform the 
CosmoMC package and the modified CAMB with the 
cosmological data including the CMB data from Planck 
with both low-/ (/ < 50) and high-/ (/ > 50) parts and 
WMAP with only the low-/ one as well as the BAO 
data from BOSS DRll and the SNIa data from SNLS. 





Fig. 4 Marginalized probability for the inflation power 
parameter q in (a) Starobinsky and (b) exponential grav¬ 
ity models, where the long-dashed, solid and dashed lines 
correspond to A = 50, 60 and 70, respectively. 


We modify two parts in the CAMB and CosmoMC 
packages. Firstly, we study the inflation model 

by adapting a standard running model to parametrize 
the primordial scalar and tensor perturbations. 


lniP«(fc) 


\nVt{k) 


In + {ris — 1) In 


+ Q!s 



In At -b rit In 




(37) 

(38) 


where kg is the pivot scale, is the scalar (tensor) 
amplitude at the pivot scale with At = rAg. Secondly, 
we examine the scalar evolution b y following the sim i- 
lar process of the MGCAMB code ( Hoiiati et al. 2011 ). 
in which Eqs. m and (Ea) are used to incorporate 
the F{R) effect. Moreover, the tensor modihcation in 
Eq. (1551) is also considered in our program. 

In Fig. m we depict the ID marginalized probabil¬ 
ity plot for the inflation power q with (a) Starobinsky 
and (b) exponential gravity models. As estimated in 
Sec. 12.11 the 2 —a allowed region indicates q ^ 0(10“^). 
In particular, we have q < 2.66 x 10“^ and 2.17 x 10“^ 
in Starobinsky and exponential F{R) models, respec¬ 
tively. 


Fig. 5 Contour plots in the planes (a) A ^ — ^Ich^ for 
Starobinsky and (b) /3“^ — for exponential gravity 

models, where the inner and outer curves represent 1 — a 
and 2 —a confidence levels, while the long-dashed, solid and 
dashed lines correspond to A = 50, 60 and 70, respectively. 


It is well-known that the F(R) gravity exists a 
ACDM limit at A (/3) —f oo with Ai?c (/3f?c) being equal 
to the cosmological constant A. We examine all the al¬ 
lowed parameter space when we fit the model parame¬ 
ters with the prior 0 < (/3~^) < 1. The constraints 

on the model parameter A“^ (/3~^) versus the physical 
density of cold dark matter (CDM) fich^ are shown in 
Fig. [5^ (b), where the long-dashed, solid and dashed 
lines represent the e-foldings of A = 50, 60 and 70, re¬ 
spectively. We note that the best-fit of the Starobinsky 
(exponential) gravity model parameter locates around 
A~^ ~ 0.75 (/3~^ ~ 0.87) with similar to the range 
in the ACDM model. These results indicate that the 
evolution history in the Starobinsky model prefers a 
significant deviation from that in the standard cosmol¬ 
ogy, so that the ACDM limit with A~^ —>■ 0 has been 
excluded within 2 —cr confidence level. On the contrary, 
the exponential type model be haves a smaller de viation 
from the standard cosmology (jCeng et al.lPOldl) . 




Fig. 6 Contour plots in the planes Emi, — flch^ for (a) 
Starobinsky and (b) exponential gravity models, where the 
inner and outer curves represent 1 — cr and 2 — cr confidence 
levels, while the long-dashed, solid and dashed lines corre¬ 
spond to A = 50, 60 and 70, respectively. 
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Table 1 Allowed regions of the baryon density Q,bh^, CDM density Q.ch'^, neutrino mass sum Errii^, inflation power 
parameter g, spectral index Ua, tensor-to-scalar ratio r, and cts with 95% CL and ^ = 1 — (1 — fi~^) in the Starobinsky 

(exponential) model with 68% CL, where the difference of the best fit between f{R) and ACDM is defined by Ax^ = 
2 2 
X/(H) ~ Xacdm- 


Parameters 

Starobinsky 

Exponential 

ACDM 

N=50 

N=60 

N=70 

N=50 

N=60 

N=70 

lOOfV? 

2.24 ±0.05 


2.25 ±0.05 

n r)q-|-U.UD 
^■^^-0.04 

2 24+IJ-U4 

^•^^-0.05 

cy QO+U.Uf 

^•^^-0.03 

9 9(4+U.Ub 

^■^U-0.03 

nj? 

nil Q + U.UU2 

U.il8 o.o5 

0.118+);-- 

0.116+^-- 

0 1 1Q+^-UUi 

u.iiy o.o5 

nil Q+O-OUi 

U.li5 Q 004 

0 11 7+^-^^-'- 

-0.004 

0.118 ±0.003 

T,m,y/eV 

< 0.302 

0-063+;;:^^^ 

“0T09l^ 

< 0.239 

< 0.264 

0 .0861;;;-^ 

< 0.211 

lOOq 

< 2.66 

< 1.86 

< 1.37 

< 2.17 

< 1.43 

< 1.10 

~ 

ris 

u.y/z o 010 

u.y/u o.oo2 

u.y/o o.ooi 

u.yuo o,ooi 

u.y/u o.oo2 

u.y/^o.oo2 

u.yuo o,oo9 

1(P7 

< 11.1 

< 6.69 

< 4.46 

< 9.42 

< 5.63 

< 3.93 

< 125 

0-8 

1 idn+u.uiiu 

l.i4U o.i39 

1 -1 qi +U.UcSS 
1.101 0.133 

1111 

0.962^;;-- 

0.963+;;-- 

0.940+;;;-^ 

A oqq+u.u:^4 

U.099 o.o59 


u.zoy o.289 

n 9Q7+u-iii 

U.ZO/Q 287 

U.ZOO o.258 

n -| 97+U.848 

-0.127 

n 1 q9+U.c}.ii 
U.10Z Qj^32 

“0T97f^ 

- 

Ax" 

-2.65 

-2.61 

-1.74 

-1.32 

-0.88 

0.26 

- 


Since the massive neutrinos play an important role 
in not only the cosmological evolutions but also particle 
physics experiments, we present the total neutrino mass 
Emjy versus the dark matter density in Fig. [51 In 
the Starobinsky (exponential) model, the allowed neu¬ 
trino mass is given by T,m^ < 0.35 (0.30) eV. Com¬ 
paring to the ACDM result of < 0.211 eV, the 
F{R) gravity results in a significant enhancement on 
the neutrino masses, especially the Starobinsky model. 
Our results for the combined f{R) scenarios are sum¬ 
marized in Table. [1] Due to the difference, listed in 
Table. [H we observe that the viable f{R) gravity mod¬ 
els in most of the examples perfectly fit the cosmic data 
(x^ < 0), especially those with the smaller e-folding N. 

4 Conclusions 

We have investigated the combined /(i?) models to 
unify inflation and dark energy scenarios. The inflation 
model is based on the generalized Starobinsky f{R) in¬ 
flation theory, R^~‘^, while the late-time dark energy 
is described by the viable f{R) gravity, including the 
Starobinsky and exponential viable f{R) gravity, be¬ 
longing to the power-law and exponential types of mod¬ 
ified gravity models, respectively. In order to exam¬ 
ine the consistency of the slow-roll parameters in the 
and the combined f{R) models, we nu¬ 
merically plot the potentials in both models and find 
that there is no difference between these two models in 
the inflationary stage but not in the post inflationary 
epoch. The difference might cause a change in the re¬ 
heating history. In addition, the behavior of the curva¬ 
ture in the Jordan frame is also changed: the curvature 
in the R^ inflation oscillates around R = 0. However, 


when we include the viable f{R) dark energy, the cur¬ 
vature keeps positive-definite and is forbidden to cross 
the critical curvature Rcr- 

The constraints from the observational data is per¬ 
fectly fitted by the combined model (x^ < 0) andpoint 
out that although the inflation can loosen the 

tensor-to-scalar ratio bound in the R^ inflation, r is still 
restricted to a tiny value by the accurate observations, 
resulting in g < 3 x 10“^, which is hardly to be distin¬ 
guished from the standard R^ inflation. We have also 
examined the model parameter of A (/3) of the viable 
Starobinsky (exponential) f{R) gravity to test the al¬ 
lowed parameter space. Although there exists a ACDM 
limit at A“^ (/3~^) 0, the Starobinsky model prefers 

^ 0.7, so that the ACDM limit is excluded within 
2 — cr confidence level, while the best-fit value in the ex¬ 
ponential gravity locates at ^ 0.85, which deviates 
from the ACDM model within 1 — cr confidence level. 
We have also shown the deviations of the neutrino mass 
sum YiTrii, between f{R) gravity and the ACDM model. 
Explicitly, we have found that the allowed could 
be released to 0.36 and 0.30 eV in Starobinsky and ex¬ 
ponential F{R) dark energy models, respectively. 
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